Spinless fermions and charged stripes at the strong-couphng hmit 
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Spinless fermions on a lattice with nearest-neighbor repulsion serve as a toy version Hubbard model, 
and have a symmetry-broken even/odd superlattice at half-filling. At infinite repulsion, doped holes 
form charged stripes which are antiphase walls (as noted by Mila in f994). Exact-diagonalization 
data for systems up to 36 sites around 1/4 filling, and also for one or two holes added to a stripe of 
length up to 12, indicate stability of the stripe-array state against phase separation. In the boson 
version of the model, the same behavior can be stabilized by addition of a four-fermi term. 
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Forty years of study have not yet produced a com- 
plete understanding of the phase diagram of the Hub- 
bard model, the simplest, nontrivial paradigm of inter- 
acting- spinfuU fermions .u The spinless lattice fermion 
modefl is a simpler and more tractable analog which re- 
tains many Hubbard-model properties, much as the Ising 
model stands in for the n-coniponent magnet in critical 
phenomena: understanding of the spinless model may 
provide fresh viewpoints of the Hubbard model, or new 
tests of known methods. Spinless models also arise nat- 
urally for ferromagnetic materials in which one of the 
spin-split bands is, completely full or completely .empty, 
such as magnetitcEl or "half-metallic" manganitesQ. 

Our aim is to promote the systematic study of this 
model's pii^e diagram, which is almost untouched in the 
liter atureOij. As a beginning, this paper argues that, in 
the strong-coupling limit, the spinless model possesses a 
phase with the quantum-fluctuating, hole-rich antiphase 
domain walls known as "stripes." Such stripes are an 
active topic in the Hubbard or t — J modelsjjitj particu- 
larly since stripes were observed in cupratescB and seem 
related to incommensurate correlations found in high- 
temperature superconductors. 

Let us take a square lattice model with Hamiltonian 
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Here Ci^ and Ci are creation/annihilation operators on site 
i, hi = c^iCi, and "(...)" counts each nearest-neighbor 
pair once. In most places we miJl consider hard-core 
bosons in parallel with fermionstZI. From here on, we 
take V/\t\ = oo so neighboring particles are simply for- 
bidden, and t is the only energy scale. (This constraint 
amounts to adopting a hardcore radius just over 1 lattice 
constant.) 

Phase diagram as function of n - Consider first the 
dilute limit, n^0.15. When V oo, the Hartree-Fock 
approximation gives absurd results; in reality the renor- 
malized interaction of two particles is of order t and a 
Bose or Fermi liquid is expected, as when the on-site re- 
pulsion f/ — > oo in the Hubbard modelt3. 

At the other extreme, the dense limit (n = 1/2, half- 
filling) admits only the two microstates with the \/2 x 



a/2 checkerboard pattern, called the "CDW" ("charge- 
density wave" ) order An Ising symmetry breaking be- 
tween even/odd lattices is exhibited, the spinless model's 
cartoon of the Heisenberg antiferromagnetic order at 
half-filling in the large- 1/ Hubbard model. 

Stripes in a hard-core model - Now consider light hole 
doping, 1/2 — n <^ 1. An isolated hole is immobile and 
gains noi-hopping energy in the CDW state (see Fig. |^). 
As MilaEj observed, a droplet including >3 holes can 
fluctuate but is still confined to a circumscribed rectan- 
gle with edges along the 45° directions, since a particle 
is prevented from hopping away from a CDW domain 
surface oriented along {11} (Fig. |l|). 

The natural way to dope holes is a "stripe" , an 
antiphase domain wall with charge 1/2 hole per unit 
length. This permits hop|S-( arrows in Fig. |l|) which imple- 
ment stripe fluctuations.llj A single stripe's path can be 
parametrized as a unique function y{x) (hopping never 
generates overhangs). Then y(x -I- 1) — y{x) — ±1, and 
the steps up or down can be represented by a string of 
corresponding "f and J, arrows. A particle hop has the 
effect 14 — > It and so the Hamiltonian jjf a single stripe 
maps exactly to the spin- 1/2 XX chainllj with exchange 
J± — —t for the X and Y spin components. That model 
is exactly soluble by a well-known mapping, whereby the 
up (down) spins map to noninteracting spinless fermions 
(empty sites), respectively, in one dimension. It follows 
that the energy of the (coarse-grained) stripe is 



dx [cto 
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where cto ~ — (2/7r)t and the stripe stiffness K — (7r/2)t. 
As usual, the sound velocity of the stripe's capillary 
waves is the Fermi velocity v = 2t/h of the ID noninter- 
acting fermions. Knowing K and v allows us to compute 
the fluctuations of the Fourier mode at each wavevec- 
tor q, as for any harmonic string: = hv/2K\q\. 

The general result is that one such "Gaussian" quantum- 
fluctuating stripe has divergent fluctuations. 



{\y{x) - 2/(0)p) = {v/2ttK) In + const. 



where v/2ttK = 2/h'^ in this case. 
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The thermodynamic, phase at 1/2 — n<C 1 could then 
be an array of stripesE3 all parallel (on average) to either 
the X 01 y axis. They have only a contact interaction, so 
the array's long-range order depends on stripe collisions, 
which surely exist since isolated stripes have divergent 
fluctuations (eq. (|l))-n 

Prior analytic workEl suggested that spinless fermions 
in d = 2, when doped away from half- filling, develop 
an incommensurate ordering wavevector slightly off from 
(tt, tt)). Expanding around the d = oo limit,Ll small dop- 
ing led to coexistence between the half-filled CDW and 
a slightly incommensurate state (but not at V/t — ^ oo). 
We conjecture these incommensurate phases, in d = 2, 
consist of stripe arrays. In the hard-core boson model 
near half- filling, in a regime < V/t < oo, the unifomi 
CDW phase is asserted to phase-separate upon dopingcil. 
The dense coexisting state is just as plausible, a priori, 
to be a strinp. array as the phase-separated state that was 
assumedtltJ. 

Stability estimates - The key question is whether (or 
when) the stripe-array is stable, compared to a phase- 
separated state in which the CDW state and the dilute 
(=hole-rich) liquid coexist. In the case of the Hubbard 
model, it was asaued that doping invariably kads to 
phase separationEa except when it is "frustrated"la by the 
long-range Coulomb force. Contrarily, it was argued that 
holes in fluctuating stripes may gain more kinetic energy 
than they would in a phase-separated stataJ. 

To decide the issue of coexistence, one first plots the 
energy per site i?'"^(n) and E^^{n) for the low-density liq- 
uid and the stripe-array, respectively, which should look 
like Fig. || for either fermions or hardcore bosons. We 
have 

E^'^^in) = i-4t)n + A^n^ + A^n^ + . . . (4) 

where the leading coefficient is the bottom of the single- 
particle bandEj With increasing density, the energy 
_B'"i(n) turns upwards and becomes small around n = 0.3 
as the hopping becomes "jammed" (neighbor sites be- 
come forbidden due to other adjacent particles). The 
matrix elements contribute with the same sign in the 
boson ground state but can't in the fermion case, so 

fermions ^ "^bosons- 

On the other hand, in the conjectured stripe array near 
half filling, 

£;^^(n)«[(7o-f 0(d)](l-2n) (5) 

where cto is the energy per unit x-length from (^; the 
second factor is the length of stripe per unit area. The 
mean stripe separation is d = (1 — 2n)~^, and (l>(d) 
parametrizes the energyjjajst per unit length from colli- 
sions of adjacent stripes, Ej so </)(d) ^ as c? — > oo. Thus, 
the chemical potential in the limit of separated stripes 
is [i* = dE^''{n)/dn\n=i/2 = (4/^^)* = 1.273i. We em- 
phasize that the form and the leading coefficients in (^) 
and (^ are the same for fermions and hardcore bosons. 
[Indeed, the fermion and boson models are identical in 



the single-stripe sector, since the V = oo constraint pre- 
vents any permutations; this identity extends to the sin- 
gle stripe with one extra hole, in which case only even 
permutations are accessiblec3.] 

There are three necessary conditions for the stability 
of the stripe array: we have strong evidence for each of 
them, frnni exact diagonalizations of systems with 20 to 
72 sites.Ea These are far too small for direct observation 
of a fluctuating stripe array, or of the coexistence of the 
liquid and dense phases; yet they are large enough to 
yield some of the parameters which the phase diagram 
can be calculated from. From here on we use units t = 1. 

The first stability condition is that stripes repel, i.e. 
(j){d) > 0; stripe attraction would suggest instability to a 
domain of liquid phase, which is scarcely distinguishable 
from a bundle of self-bound stripes. In both the fermion 
and boson cases we diagonalized L x L' systems, for L = 
4, and L' — 6, 8, 10, as well as L = L' = 6, doped with 
L holes so that two stripes run in the short direction, 
and thus d = L' /2. Define a stripe interaction per unit 
length 0off(^72) = {E - 2E''^'{L))/2L, where E^'^'iL) is 
the energy of one isolated stripe of finite (even) length 
L. {E^*^'^{L) is calculated exactly using the ID spinless- 
fermion representation.) Indeed, 4>csid) was positive and 
(for L = 4) decreasing with d. j— . 

Next, a stripe can contain extra holesjlj which move 
as quasiparticles with an excitation gap A. The second 
stability condition is 

A > -/I* (6) 

If not, further doping would add holes to existing stripes 
rather than form new ones, again suggesting a ten- 
dency to form phase-separated droplets. From diago- 
nalizations we measured the excitation energy of one 
added hole,EZl A(L,L'), with even L and odd L' in 
the range [4, 8] for fermions, or [4, 12] for bosons; as 
mentioned above, this is strictly independent of statis- 
ticsE3. From this we extrapolated, first L' — )■ oo using 
A(L, L') = A(L) + AA(L)e~'^'/'(^) and then L oo us- 
ing A(L) = A + Sa/L. We found A = -0.65(5), which 
comfortably satisfies (^. 

We also analyzed the two-hole energy A2(L), for L — 
4, 6, 8, 10 only; with two holes, the L dependence is more 
like 1/2-^ than 1/L. Extrapolating to L = oo, yielded 
A2 = —1.42(7) for bosons and —1.44(9) for fermions. 
Note A2 — 2A « —0.1, i.e. hole binding is insignificant 
when L < 10; we think it is a real effect in a large system, 
since holes on a stripe can be collected into a (hole-free) 
vertical segment of the stripe. (Since the stripe's 90° 
kinks cost energy, an array of parallel stripes will still be 
the thermodynamic phase in a large system.) 

The third condition is the crucial one: as shown in 
Fig. 1^, the dense phase coexisting with the liquid must 
not be the CDW, which would preempt a stripe-array 
phase. That is, 

// > (7) 
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where /z is the slope of the trial liquid-CDW tie-line 
tangent to and passing through (n = 1/2, E = 0). 

To test (0), the equation of state is required. 

We exactly diagonalized all rectangular lattices with 
L,L' > 4 and LL' — 20 to 36 sites, and with occupa- 
tion in the range 0.20 < n < 0.3, and fitted the results 
to (|). We obtained (^2,^3) ~ (9.45 ± 0.6,5 ± 2) for 
bosons and (11.25±0.6, — 1±2) for fermions. This implies 
Ai^^'bo.ons = 1-33(2) and M^^fcrmions = 1-25(2). (Here the 
errors are estimated by varying the subset of data used 
for the fit.) For either bosons or fermions, coexistence 
with the CDW would occur at ric ~ 0.24. 

Hence, stripes are unstable in the boson case and (very 
likely) stable in the fermion case, but close enough to the 
boundary in either case that the balance can be tipped 
either way by the small perturbation tc, discussed later. 

Exotic states? - Like the Hubbard model, the spinless 
fermion model may be extended by adding other hop- 
ping terms to the Hamiltonian, which might stabilize ad- 
ditional phases. Many of these terms have the form 

- ix{c^ jCi + c'tCj)hk, (8) 

where {i,j,k) are three different sites arranged as in 
Fig. H, and x is "a" , "6" , or "c" for the hops shown in the 
corresponding parts of Fig. ^. For example, when V is 
large but finite, hops are possible to a neighbor's neighbor 
with ia = ib = —t^jV as in Fig. ||(a) and (b), analogous 
to similar terms of order jU when the t-J model is de- 
rived from the Hubbard model. This spinless analog of 
the t-J model, in which virtual states with neighbor pairs 
are projected out, will be the natural starting point to 
study phenomena at large (but not infinite) V ^ e.g. the 
mobility of lone holes. 

In the fermion model, one could artificially take tb ^ t 
(analogous to J ^ i in the t-J model); then the term 
Fig. [s] (b) naturally favors superconductivity. Namely, 
fermions form tightly bound p-wave pairs, separated by 
V^; these composite bosons hop with bandwidth 8tb, and 
Bose-condense in the usual fashion. Thus tb is analogous 
to negative U in the Hubbard model, in that supercon- 
ductivity is put in "by hand" . But it is a plausible spec- 
ulation that, in the highly correlated liquid at n sa 0.25, 
BCS superconductivity appears even with 4 ^ t. 

Finally, consider the hopping ic of Fig. || (c) which 
just modifies the amplitude of already possible hops. 
This tends to stabilize (destabihze) stripes according to 
whether ic has the same (opposite) sign as t, since every 
allowed hop in a stripe is surrounded by particles on all 
four possible "fc" sites. (Compare Fig. g(c) with Fig. H). 
Hence the stripe energy (Tq and /i* get multiplied by a fac- 
tor (1 -f 4tc/t). On the other hand, assuming that each 
"fc" site is about 1/3 occupied in the liquid at n = 0.25, 
it follows that i?''i(n) is multiphed by about (1 + 4^/3^). 
If so, the critical perturbation where — /i* (so the 
stripe phase appears or vanishes) is only tc/t « 0.02 for 
bosons or —0.007 for fermions, using our values of fj}"'^ 
quoted above. 



Discussion - To establish the occurrence of stripes in 
the V — oo system, we addressed, by exact diagonaliza- 
tions, (i) stripe-stripe interactions; (ii) the energy of a 
single hole, as well as hole-hole interactions, on a stripe; 
and (iii) the medium density liquid regime, (iv) hole- 
hole interactions on a stripe The enormous reduction of 
Hilbert space due to the nearest-neighbor exclusion (at 
V = oo), as well as the lack of spin, permits numerical 
explorations at system sizes much larger than would be 
possible in the Hubbard model - vital not just for study- 
ing stripes, but any microscopically inhomogenous states. 
Monte Carlo simulation of the strip|e-| phase is straight- 
forward for the hardcore boson case.E3 Boson results are 
valid for fermions too, when the stripe separation d is 
large and the density of "extra holes" on each-stripe is 
low, since particles do not exchange in this limited. As for 
the fermion case, the new "meron-cluster" Monte Carlo 
algorithm cancels the sign problem for a limited class 
of models including spinless fermions, but not the basic 
Hubbard modelEa. 

The ultimate aim of microscopic simulations should be 
to extract macroscopic parameters, e.g. the stripe stiff- 
ness K or the stripe contact repulsion. This is more 
straightforward than in the spin- full (Hubbard or t-J) 
case, where the inter-stripe domains contain gapless spin- 
wave excitationstj. These parameters may be input to 
analytic explorations of the interesting anisotrooic con- 
ductivity of the quantum-fiuctuating stripe arraycj, also 
simpler in the spinless case. 

More broadly, it is a challenge to test for the exotic 
phases we mentioned in connection with Fig. ^. In the 
medium-density regime n w 0.2, strong correlations of 
some sort are essential to minimize the hopping energy. 
These might be prosaic, e.g. a -s/S x ^/b superlattice, 
but the following possibilities are realizable, in principle, 
even in a spinless model: (i) orbital magnetism (sponta- 
neous circulating currents around plaquettes) ; (ii) p-wave 
superconductivity (see the speculations on Fig. ^(b)); 
or (iii) the analog of spin-charge separation, the spinon 
being replaced by a spinless particle that carries Fermi 
statisticfS. If it transpires that such states are hard to 
stabilize without spin, that would shed additional light 
on the Hubbard model; contrariwise, if they are stabi- 
lized, they may be easier to study in the spinless case, 
free from any background of low-energy spin excitations. 

A crude comparison may be made of the spinless- 
fermion model with the infinite-?/ Hubbard model in the 
dilute regime. In that case, each fermion excludes one 
site (its own) from half oi the other fermions, not count- 
ing the exclusion built in by Fermi statistics. In the 
present spinless model each fermion excludes four sites 
from all other fermions, so in a sense the hole-rich metal 
phase is "jammed" eight times more effectively than in 
the Hubbard case. We expect, then, that kinetic-energy- 
stabilized stripes are far more robust in the present model 
than in the large-[/ Hubbard case. In fact they are prac- 
tically marginal in the present model, so this encourages 
the opinion that stripes are not stable in the short-range 
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Hubbard (or t~J) model. 
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stripe lone hole hole on stripe 4^" interface 
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FIG. 1. Hard-core particles (filled circles) must form a 
checkerboard at half-filling in the V = 00 limit. No hopping is 
possible by the lone hole ( x mark) or 45° interface (dot-dash 
line), but hopping is easy (arrows) along the stripe. The 
stripe's edges are shown by shaded lines and its spin-chain 
representation is below. One extra hole is shown on the stripe. 
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FIG. 2. Schematic of energy E{n) as a function of filling n, 
assuming < fj* . Dotted segments are metastable states. 
Thin lines are drawn with slopes /i* and /x'"^. 
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FIG. 3. More hopping terms allowed in the V — ca Hilbert 
space, from site i to site j when site k is occupied (Eq. 
These terms might modify the phase diagram (see text). 
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